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Abstract. We compare the properties of subwavelength imaging in the visible 
wavelength range for metal-dielectric multilayers operating in the canalization and 
the resonant tunnelling regimes. The analysis is based on the transfer matrix method 
and time domain simulations. We show that Point Spread Functions for the first 
two resonances in the canalization regime are approximately Gaussian in shape. 
Material losses suppress transmission for higher resonances, regularise the PSF but 
do not compromise the resolution. In the resonant tunnelling regime, the MTF 
may dramatically vary in their phase dependence. Resulting PSF may have a sub- 
wavelength thickness as well as may be broad with multiple maxima and a rapid phase 
modulation. We show that the width of PSF may be reduced by further propagation 
in free space, and we provide arguments to explain this surprising observation. 
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Figure 1. a. Schematic view of the multilayer consisting of N periods of thickness 
A. Altogether the multilayer contains 27V + 1 layers of silver, or dielectric x\ b. 
Schematic view of an experimental set-up for the characterisation of imaging through 
the multilayer [7]. 



1. Introduction 

In 2000, Pendry pQ has shown that a 40nm silver slab is capable of near filed imaging 
with resolution beyond the diffraction limit in the near UV wavelength range. The 
operation of a similar but asymmetric [2] flat lens has been later experimentally 
verified [3j Hj. In these works, the mechanism leading to superresolution is related 
to the appearance of surface plasmon polariton (SPP) mode, which enables the 
transfer of evanescent part of the spatial spectrum through the slab. Other ways of 
achieving superlensing include the use of photonic bandgap materials showing negative 
refraction [5]. The negative refraction is similar as in left-handed metamaterials 
discussed forty years ago by Veselago [6] who studied theoretically the properties of 
materials with a negative real part of both permittivity and permeability. 

Moreover, also the evanescent waves are subject to amplitude enhancement in 
the left-handed materials. Therefore, the imaging through a slab with a Veselago 
material may be, at least theoretically, truly lossless and aberrationless. For the 
Pendry's superlens, both losses as well as the cut-off wavelength of the SPP mode 
limit the operation range to subwavelength distances only [8J [9]. This obstacle has 
been tackled to some degree in several ways. One of them is by breaking the silver 
slab into a metallic multilayer [2j [10] with effective anisotropic properties and strong 
coupling of the SPP modes between neighbouring layers. Otherwise the loss may 
be compensated by the use of media with gain [HI [12]. Unfortunately, SPP modes 
introduce singularities into the modulation transfer function (MTF) [131 [H], resulting 
in deteriorated imaging j2j [151 EE 02] • 

The aim of this paper is to characterise the sub-wavelength near-field imaging 
properties of superlenses operating in two other, recently discussed regimes. These 
operation regimes are the resonant tunnelling [THl HH] and canalization [2Ql EH [22] . 
Unlike the previously discussed regimes, which used SPP or negative refraction, resonant 
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tunnelling and canalization rather depend on the band structure of the metal-dielectric 
one-dimensional photonic crystal, its effective properties, and the coupling properties of 
the external layers. We analyse structures with a high transmission, with thicknesses 
ranging up to several wavelengths - working under resonant tunnelling, or thinner than 
wavelength but with very regular responses - working under canalization. We also 
address the role of losses and finite layer thickness for the canalization regime. In our 
analysis we characterise these superlenses with the point spread function (PSF) [T3l[T4] . 
and their FWHM and standard deviations. The potential applications of these elements 
are similar to other superlenses and include eg. nanolithography [4j, and proximity 
printing [TT] . 

2. Background 

2.1. Numerical method 

Our numerical analysis combines two elements: 1. the transfer matrix method 
(TMM)[23j, which we use to examine the transmission of a single plane- wave through 
the multilayer, and 2. the theory of linear shift-invariant systems (LSI) p~3j[H], which 
provides the framework for the description of the multilayer with the modulation transfer 
function (MTF) or the point spread functions (PSF). In this section we review the 
principles of both elements. Appart from TMM, we also include finite difference time 
domain (FDTD) [24] simulations which are performed using a freely available software 
package MEEP with subpixel smoothing for increased accuracy [25j . 

Linear shift invariant scalar systems (LSI) are well known in different fields of 
Optics. Here we apply this model for the description of optical multilayers in a situation 
when they act as imaging elements for coherent monochromatic light. 

Let us introduce some denotations that are used throughout the entire text. 
Figure [jji shows a schematic view of a metal-dielectric periodic multilayer. The structure 
is surrounded by air, and consists of periodically arranged metallic (Ag) and dielectric 
layers (made of a dielectric denoted as x). The structure has the total thickness 
L, contains TV periods with the pitch A, and single layer thicknesses of d,A g and d x , 
respectively. First and last layers have the thickness dependent on parameter p. The 
initial layer of x has the thickness of p(A — cU^)? an d the last has the thickness of 
(1 — p)(A — dA g )- The incident light is coherent and monochromatic with the free-space 
wavenumber k = 2tt/X. The time dependence of the complex fields is exp(—ik ct). 

In-plane imaging through a layered structure consisting of uniform and isotropic 
materials represents a LSI, for either TE or TM polarisations. Linearity of the system is 
the consequence of linearity of materials and validity of the superposition principle for 
the optical fields. Shift invariance results from the assumed infinite perpendicular size 
of the multilayer and the freedom in the choice of an optical axis. A scalar description 
is valid for the TE and TM polarisations in 2D since all other field components may be 
derived from E y or H y) respectively. For instance, for the TM polarisation, the magnetic 
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field H y {x, z) may be represented with its spatial spectrum H y (k X) z) 



Hy{X ) Z) 



+00 



H y {k X ) z)exp(ik x x)dk X) 



(i) 



where, at least for lossless materials, the spatial spectrum is clearly separated into the 
propagating part k 2 x < k^e and evanescent part k 2 x > k^e. 

LSI can be characterised either using the modulation transfer function (MTF), or 
equivalently with the point spread function (PSF) [TSKH]. The MTF is simply the ratio 
of the spatial spectra of the output and incident fields and therefore corresponds to the 
amplitude transmission coefficient of the multilayer 



H y (k x , z = L) = MTF(k x ) • m nc (k xi z = 0). 



(2) 



We note that due to reflections, the incident field Hy(k X) z — 0) differs from the total 
field H y {k X) z = 0). 

MTF is obtained with the TMM method. For this purpose we analyse separately the 
transmission of plane- waves with different real values of the wavevector component k X) 
corresponding to both propagating and evanescent waves. In each layer, the wavevector 



component k x is conserved, and k l z = iy/c^ — k 2 . The magnetic field in layer / may 
be written as 

H l y (x, z) = A l + exp(i(k x x + k[z)) + A l _exp(i(k x x - k l z z)), (3) 

where the amplitudes A l + and A l _ are matched at layer boundaries to satisfy the 
continuity of both H y and E x . This leads to the transfer matrix relation between 
subsequent layers 
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where z\ is the position of the boundary, and r is the Fresnel reflection coefficient for H y 
at the boundary. Solving the linear system of equations from eq. Q gives the complex 
transmission and reflection coefficients of the multilayer, and together with eq. ^ allows 
to reproduce the field inside the structure. We will publish separately the details of the 
numerical method which fall beyond the scope of this paper. Let us only note that the 
1. TMM further simplifies for periodic structures, 2. using a moving coordinate system 
for the amplitudes A l _ and A l + improves stability of the TMM for evanescent waves, 3. 
a full 3D analysis is more involved [26] and falls beyond the scalar model, 4. the present 
description is compatible with plane- wave waveguide solvers, eg.[27j, which may easily 
serve as the source field generator. 

PSF is the inverse Fourier Transform of the MTF and has the interpretation of 
the response of the system to a point signal 5(x). The response to an arbitrary input 
Hy nc (x, z — 0) can be further expressed as its convolution with the PSF 

H y (x, z = L) = H^ nc (x, z = 0) * PSF(x) (5) 

While the MTF and PSF represent the imaging system in an equivalent way, and 
MTF is easier to calculate, there exist obvious advantages of using PSF. PSF can be more 
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Figure 2. Real and imaginary part of the permittivity of silver [28 , and the 
corresponding skin depth. 

straightforwardly interpreted, as it provides visible information about the resolution, 
loss or enhancement of contrast, as well as the character of image distortions. We 
use the standard deviation and full-width-at-half-maximum (FWHM) to compare the 
characteristics of PSF. However, we have to note that only for Gaussian PSF and input, 
the output given from eq. ^ has the width (variance) equal to the sum of the widths 
of PSF and input. Normally, on this basis it is possible to link the resolution with the 
thickness of PSF. In the present analysis we also deal with PSFs very different from a 
Gaussian function. 

Our present analysis is limited to transmission through multilayers surrounded by 
air and does not include a particular source and detector. An experimental verification 
of these results is possible in a set-up presented recently by Fabre at.al [7J for the 
characterisation of superlenses and schematically shown in fig. [TJd. A realistic element 
for use in eg. proximity lithography or for coupling between a plasmonic circuit and 
an optical fiber requires certain changes to our current assumptions. It would include 
different materials on both sides of multilayer, and have the input and output planes 
placed further from the structure. Its design would be preferably based on numerical 
optimisation of the structure. In this paper instead, we have deliberately chosen to rely 
on theoretical models represented by the two clearly defined regimes. 

2.2. Review of the resonant tunnelling and canalization regimes 

Resonant tunnelling and canalization are the two regimes compared in this paper. 
Both of them enable a metal-dielectric multilayer to perform the role of a near-field 
subwavelength imaging device operating at visible wavelengths. 

Resonant tunnelling has been proposed and demonstrated already in 1998 by 
Scalora, Bloemer et. al PJ2 [29], who have shown that in the pass-band of the metal- 
dielectric stacks around the visible wavelengths, light penetration into metal is not 
limited to the skin depth and can reach distances of the wavelength order. This 
phenomenon called transparency of metals has been recently recalled [18] in the context 
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Figure 3. MTF and the reflection coefficient for the multilayer with e x = 1 — Re(cAg), 



Re(e Ag )/e 



-dAg/d x , and p = 0. Resonant condition for canalization is achieved for 



2L/X G Z. a. lossless model with a large number of periods N — 2000, €A g = —12.88; 
b. lossless model with a realistic number of periods TV = 20, €A g — —12.88; c. lossy 
model with a realistic number of periods N = 20, €A g — —12.88 + 0.48z. 



of near-field subwavelength imaging. Some other interesting features of this regime 
have been also indicated [18], such as the possibility to tailor the transmission band 
in a wide range, or the possibility of enlarging the structure with the surprising effect 
of increasing the transmission. It has also been shown that the transmission depends 
strongly on the first and last layers which couple and decouple light to and from the 
periodic structure [18J. 

The second regime to be included in our comparison, named canalization has been 
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introduced Belov et. al. [21J and experimentally demonstated by Ikonen et al. [22] 
for a capacitively loaded wire medium. In 2006 Belov and Hao [20]j extended the 
concept of canalization to metal-dielectric layered structures. We will only refer to 
such type of structures in the present paper. The canalization regime [20] enables 
subwavelength imaging at visible frequencies through metal-dielectric multilayers. The 
multilayer consists of layers sufficiently narrow to be properly characterised with the 
homogenisation predicted by the effective medium theory (EMT). In the canalization 
regime, the multilayer needs to have the effective permittivity tensor of the anisotropic 
form with e± = oo, and e\\ — 1, and at the same time to fullfill the Fabry- Perot condition 
for minimal reflections 2L/A G Z. This leads to simple geometrical conditions of 
£Ag/t x = —dA g /d x and e x = 1 — €a 9 - To fulfil the last condition, one has to neglect 
the losses in silver, and assume that e^, and e x are real but have a different sign. In 
fig. [2] we show the dispersion of 6a 9 for visible and near-UV wavelengths [28]. Indeed the 
magnitude of imaginary part of 6a 9 is significantly lower than the magnitude of its real 
part, while 1 — Re(eA 9 ) falls in the range of most known dielectric materials. Therefore, 
the canalization regime may be realised with use of silver and a variety of dielectric 
materials. 

3. Canalization 

Canalization regime can only be achieved when the silver and dielectric layers have 
real permittivities, and 6A 9 /e x — —dA g /d x . Since permittivity of metals is complex, this 
relation can only be approximately fulfilled. Moreover, the validity of EMT requires that 
the layers were small with respect to wavelength A << A. In this section we investigate 
the practical importance of both assumptions. For the purpose of this analysis, let us 
assume that the wavelength is in the middle of the visible range A = 550nm, when 
the permittivity of silver is equal to 6a 9 = —12.88 + 0.48z. The dielectric x needs to 
have the permittivity of e x = 13.88, and dA g /d x « 0.928. In fig. [3] we characterise 
the MTF of such a multilayer, as well as the corresponding reflection coefficient as 
a function of L/A, and k x /k . Both amplitudes and phases of the MTF are shown. 
The values of the wavevector component parallel to the boundaries k x range from 
(normal incidence), through < k x < k (oblique incidence from air), up to larger 
values k x > k (evanescent waves) important for sub- wavelength imaging. In fig. [3] we 
compare the MTFs calculated for an unrealistic number of periods TV = 2000 and no 
losses in metal (fig. [3^i), with a more realistic case N = 20 (fig. [3)3), when the EMT 
applies only approximately. Finally, in fig. [3)3, we included losses into the calculation. 
Canalization requires that 2L/A G Z. This is the condition to minimise reflections from 
the Fabry- Perrot slab with effective anisotropic properties independently of the angle of 
incidence. Several conclusions can be drawn from the MTFs in fig. [3j The finite width 
of the layers primarily affects the high FP resonances, shifts their location towards 
smaller values of LjX and introduces a certain dependence of the resonant condition on 
k x /k . Further, the transmission bandwidth is reduced by losses (fig. |3fc). Moreover, 
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Figure 4. Transmission of a plane wave through a metal-dielectric multilayer 
operating in the canalization regime. a,d. normal incidence, k x = 0; b,e. oblique 
incidence, k x = 0.5/co; c,d. evanescent wave, k x = 1.5fco5 Parameters of the multilayer 
are the same as in fig. [3^. with L/X = 0.47 (abc) or L/X — 0.87 (def). 



losses introduce a rapid phase variation in the vicinity of k x = fc , which must result in 
deteriorated imaging. In practise, only the first two resonances allow for sub-wavelength 
imaging, therefore limiting the thickness of the slab L to subwavelength distances. 

Figs. [3^l-c indicated significant differences between the MTFs representing the 
canalization in its ideally formulated concept, and its implementation with realistic 
assumptions. Nevertheless, the distribution of field transmitted through the structure 
(see fig. [Jt-f) proves that also the realistic imaging device operates according to 
expectations. In fact, a plane wave under normal incidence (fig. [ikd), oblique incidence 



Sub-wavelength imaging in the canalization and resonant tunnelling regimes 



9 



CANALIZATION {lossless model, 1st resonance) CANALIZATION {lossless model,2nd resonance) 




Figure 5. PSF of the lossless multilayer with same parameters as in fig.[3]3. operating 
in the canalization regime a. at the L/X = 0.48 resonance; b. at the L/X = 0.88 
resonance. For comparison, the PSF for the TE-polarisation and for the diffraction in 
air at the same distance are also shown. 
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Figure 6. PSF of the lossy multilayer with the same parameters as in fig.[3j3 operating 
in the canalization regime a. at the L/X = 0.47 resonance; b. at the L/X = 0.87 
resonance. For comparison, the PSF for the TE-polarisation and for the diffraction in 
air at the same distance are also shown. 



(fig. |4pe), and an evanescent wave (fig. |4pf) are all transmitted with high efficiency. 
Moreover, the distributions of H y and E x inside the structure, all represent the same 
FP resonant condition of L/X ~ 0.5 or L/X ~ 1, which should not depend on the angle 
of incidence. For these suboptimal but realistic conditions, both for the propagating 
and evanescent waves, the multilayer is therefore characterised with a similar Bloch 
vector, and operates as a FP plate at resonant conditions independently on k x /k . On 
average, the phase of both H y and E x increases in the direction of propagation for the 
propagating waves with the same slope as in air, indicating the same effective phase 
velocity. However, locally the phase of H y propagates backwards in the silver layers. 

Resolution of the imaging system normally can be deduced from the width of the 
PSF. In fig. [5] we present the PSF corresponding to the first and second FP resonances, 
when the losses are neglected. In order to determine the actual location of the resonances 
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we minimised the ratio between reflection and transmission for normal incidence. PSF 
is calculated from the MTF in fig. [3)3., for A = 550nm. In Table [T] we compare the 
FWHM and standard deviation of the PSFs presented in this paper. The thickness of 
the multilayer is equal to L — 0.48 A = 264nm, or L = 0.88 A = 484nm, for the two 
resonances respectively. The central maximum of the \PSF\ 2 has the FWHM equal 
to 23nm or 29nm (see Table [T]), allowing for imaging with sub- wavelength resolution. 
The PSFs for free-space propagation along the same distances are much broader (see 
figs.[5^b). The PSF for the first resonance is approximately Gaussian in shape. PSF for 
the second resonance (fig. ^jp) has approximately twice lower amplitude and is distorted 
by the presence of a background structure resulting from a SPP mode that appeared 
for large k x /k in the MTF. In the presence of losses, this mode has been suppressed 
(see fig. [3)3) resulting in the regularisation of the corresponding PSF. Figure |6^b show 
the PSF for the first two FP resonances, with losses included into calculation. The 
apparent regularisation of the PSF by the presence of losses appeared at the price of 
reduced amplitude and a slower decay far from the center (this is also visible from the 
large values of standard deviation in comparison to FWHM in Tab (I). However, the 
resolution measured with FWHM is only little affected. 

We consider the approximately Gaussian shape of the PSF as a major advantage 
of the canalization regime. This highly desired feature reflects the regular shape of the 
MTF. In particular, imaging in the canalization regime does not involve plasmons, and 
for L/X = 0.5, the MTF does not experience significant enhancements for evanescent 
waves. Moreover, its phase behaves smoothly near k x /ko ~ 1. In the lossless case, phase 
of the MTF is approximately constant in the neighbourhood of that point. Losses affect 
this situation only in the very close vicinity of k x /k w 1, and a local phase modulation 
around that point has no major influence on the shape of PSF. However, due to the way 
in which losses affect the MTF, it will be difficult to utilise higher order FP resonances, 
and therefore construct multilayers thicker than the wavelength. 

4. Resonant tunnelling - broadband transparency 

We begin by stating the known differences between canalization and resonant tunnelling. 
Resonant tunnelling [19l [18], as opposed to canalization, allows for broadband 
transmission of light. The transmission window may be tuned in a large range 
and corresponds to the pass-band of the photonic crystal. Transmission has no 
straightforward dependence on homogenisation, and the layers may be thicker than 
in the canalization regime. However, the first and last layers perform a role similar to 
anti-reflection coatings and have a major influence on the transmission. They couple 
the incident wave to a Bloch wave propagating through the periodic structure, and 
decouple it in the end. No special relation between the permittivity of metal and 
dielectric components is required, leaving much space for numerical optimisation. 

We focus on the imaging properties of multilayers operating under resonant 
tunnelling. Let us consider a multilayer consisting of silver and a dielectric with 
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Figure 7. MTF of the multilayer operating in the resonant tunnelling regime described 
by the following parameters: TV = 30, p = 0.5, dAg = 20nm, e x = 8 and a. A = 55nra; 
b. A = 80nra. 



RESONANT TUNNELLING { X =450nm,A =55nm) RESONANT TUNNELLING (X =550nm,A =80ram) 
4.0 1 , , , , I 1.0, — , , , , , i 
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Figure 8. PSF of the multilayer operating in the resonant tunnelling regime a. for 
the same multilayer as in fig. [7^, at A = 450nra; a. for the same multilayer as in 
fig. [7]3, at A = 550nra; PSF is shown for the TM, and TE polarisations as well as for 
the diffraction in air at the same distance. 
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Figure 9. a. Amplitude and phase of the MTF of the multilayers operating in the 
resonant tunnelling regime from fig. [7^ (cross-section at A = 450nm) and fig.[7]3 (cross- 
section at A = 550nm); b. Phase of the central part of PSF from fig. [8]o. 
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Figure 10. Resonant tunnelling regime: ab. A = bbnm A = 450nm; cd. A = 80nm 
A = 550nra. Figures include the instantaneous distribution of H y for CW illumination 
in the steady-state, obtained using FDTD when the source is a Gaussian beam (ac), 
or a point source (bd) put bum away from the surface of the multilayer. 



an arbitrarily selected permittivity of e x = 8. In fig. [7] we present the wavelength 
dependence of the MTF and reflection coefficient for two multilayers. Both include 
N = 30 periods, with Ag layers thickness close to the skin depth d,A g = 20nm, and either 
d x = 35nm, or d x = 60nm. The losses and dispersion of silver are taken into account. 
Following the suggestion from [18], we have chosen the external layers to have the 
thickness such that p = 0.5. This value optimised the transmission in our simulations, 
as well. We emphasise that the general characteristics of the MTFs obtained with a 
wide range of different parameters are very similar to those analysed in fig. [7} These 
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Figure 11. a. PSF of the multilayer operating in the resonant tunnelling regime 
described by the same parameters as in fig. [8)3, however with the position of the 
incidence plane shifted from the boundary of the structure a. by 20nra a. by 225nm. 



MTFs are typical for the resonant tunnelling. 

The results from fig. [7^ib indicate the existence of a broad spectral region of 
transparency at 350nm — 550wn or 450nm — 760nm, respectively. Tuning to the desired 
range is based on the simple choice of thicknesses d x . Transparency windows extend 
to high evanescent waves k x /ko ~ 3, however with certain restrictions. Transmission of 
evanescent waves is accompanied by high reflections. The corresponding MTF is not 
a uniform function of k x /ko and experiences several maxima and a phase modulation. 
Finally the phase is not continuous at k x /k = 1. These features may result in a 
distorted shape of PSFs as compared to the previously discussed canalization regime. 

In fig. [8^1 and fig. [8)3 we present two selected PSFs for the resonant tunnelling regime. 
PSF form fig. [S^t corresponds to the MTF from fig. [7^i calculated for the wavelength of 
A = 450nm, while PSF form fig. [8)3 corresponds to the MTF from fig. [TJd calculated 
for the wavelength of A = 550nm. Both PSFs are significantly narrower compared to 
the corresponding PSF for propagation in air, or to the PSFs for the TE polarisation. 
While, in its centre the PSF in fig. [8^1 resembles a narrow Gaussian function, it also 
contains a slowly decaying backround. The PSF in fig. [8)3 has an interesting structure 
with two maxima and a periodic modulation. The astonishing difference between PSFs 
from fig. [8^l and fig. [8)3 results from the very different phase slopes of the corresponding 
MTFs, which we compare in fig. [9^i. The broadening of PSF effects from the large slope 
of phase in the evanescent part of the spatial spectrum. 

The two maxima of PSF from fig. [8)3 are spatially separated by around 3A. Such a 
shape of PSF could indicate the possibility of using the considered structure as a beam- 
splitter. However, this PSF is not a simple superposition of two Gaussian-like functions 
but also is subject to a rapid phase modulation (see figj^)). In effect it transforms the 
incident wavefront in a way which is different than simple splitting and broadening. 
To illustrate this, we show in fig 10 how the corresponding multilayer transmitts a 
narrow Gaussian beam comparing it to the transmission of a beam incident from a 
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RES. TUNNELLING A=55ftm ; -\=450ftm {Normal incidence k x /k =0) 
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RES. TUNNELLING A=55nm,A=450wm (Oblique incidence k x /k =0.5) 




RES. TUNNELLING A=55nm,A =450iim (Evanescent wave fe^/fc,, =1.1) 
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RES. TUNNELLING A=80nm,A =55Dnm (Oblique incidence k x /k =0.5) 




RES. TUNNELLING A =80iim, A =550nm( Evanescent wave kjha =1.1) 




Figure 12. Transmission of a plane wave through the multilayer operating in the 
resonant tunnelling regime. a,d. normal incidence, k x /ko = 0.; b,e. oblique incidence, 
kx/ko = 0.5; c,f. evanescent wave, k x /ko = 1.1; Parameters of the multilayer are the 
same as in fig. [7^. for A = 450nm (abc) or in fig. [7]o. for A = 550nm (def). 



point source. This simulation is based on the finite difference time domain method [21] 
(FDTD). FDTD Simulations were performed using a freely available MEEP software 
package with subpixel smoothing for increased accuracy [25j. The TMM is a quasi- 
analytical method, while accurate FDTD simulations of metallic layered structures are 
difficult [30] and may not reproduce the fine details of the field. Nevertheless our FDTD 
simulations come in good agreement with TMM. 

The conditions for the resonant tunnelling regime correspond to those in fig. [8^b. 
with a CW source put very close 5nm) to the surface of the multilayer. In contrast 
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Table 1. Resolution in the canalization and resonant tunnelling regimes. Emphasised 
values indicate the PSF widths below A/2 and structure thicknesses over A. 



FWHM 



Description 



Fig. 



A [nm 



\PSF\ 2 



A 



'\PSF\ 2 

A 



Canalization, lossless, 1st resonance 
Canalization, lossless, 2nd resonance 
Canalization, lossy, 1st resonance 
Canalization, lossy, 2nd resonance 



5a 

53 
6 3 



550 
550 
550 
550 



0.042 
0.066 
0.052 
0.074 



0.019 
0.049 
0.060 
0.11 



Res. tunnelling A 
Res. tunnelling A 
20nm off-plane 
225nm off-plane 



55nm 
80nm 



83 



TT 



450 
550 
550 
550 



0.17 

3.89 
3.85 
2.18 



0.99 
1.76 
1.56 
0.85 



to the point source, the Gaussian beam is neither split nor broadened (figs. 10 3d). 

This difference deserves a further explanation. The decomposition of PSF from 
fig. [8)3 into spatial harmonics differs significantly from that of a Gaussian function 
with a similar width. The Gaussian beam is composed mainly from the propagating 
spatial harmonics, while the broad shape of PSF (fig. [8)3) results from the rapid phase 
modulation in the evanescent part of the spatial spectrum (fig.|9^.). For similar reasons, 
the PSF from fig. [8)3 is also highly sensitive to the experimental conditions. For instance, 
an additional propagation in air either before or after the multilayer changes the shape 
of the PSF significantly by filtering out the evanescent part of the MTF. In fig. [TT] we 



show how the PSF is modified with distance. Already after 20nm the side-maxima are 
attenuated to the size of the central maximum (figfTT^i). After the distance of 225nm the 
PSF becomes approximately Gaussian. Its width is reduced to the central maximum 
providing a qualitative agreement with FDTD simulation for the Gaussian incidence. 

For completeness, in fig. [12] we present the internal distribution of electric and 
magnetic field for a plane wave transmitted through the multilayer. The multilayer 
operating in the resonant tunnelling regime is the same as in figs. [7)3, and [8)3. Fig. [12^, b, 
and c, show respectively the transmission of a normally incident, inclined or evanescent 
wave, through the whole structure. In contrast to fig.[4]for the canalization regime, there 
is no indication of a FP resonance this time. Indeed both transmission mechanisms are 
different, and as we show provide different PSF characteristics. 



5. Conclusion 

We compare subwavelength imaging in the canalization [20] and resonant tunnelling [T8] 
regimes. 

The canalization regime is based on the homogenisation of the multilayer and 
enables transmission for resonant frequencies or slab widths. The conditions for 
canalization can not be accurately achieved in presence of losses. We show that in 
the lossless case, the phase of MTF is approximately constant in a wide range of k x /k$. 
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Losses affect this situation mainly in the very close vicinity of k x /k$ « 1, and a local 
phase modulation around that point has no major influence on the shape of PSF. Losses 
primarily decrease the amplitude of the PSF and suppress the possibility of using higher 
FP resonances. Losses also remove the background noise in the PSF, and slightly affect 
the resolution. In effect, the PSF is approximately Gaussian in shape, which is a major 
advantage of the canalization regime. 

Resonant tunnelling enables broadband transmission through multilayers with 
thickness ranging up to several wavelengths. We demonstrate that the resulting PSFs 
may be distorted by sidelobes, and slow attenuation far from the central position. These 
effects are not entirely suppressed by losses. Such behaviour results from the shape of 
corresponding MTF, which is not a uniform function of k x /k^ ) experiences multiple 
maxima and a phase modulation, and finally the phase of MTF is not continuous at 
k x /k = 1. Nevertheless, subwavelength imaging is still possible. We demonstrate a 
situation, when the PSF becomes narrower with growing distance from the point source. 
This effect may be understood from the rapid phase modulation of both MTF and PSF. 
Further from the multilayer, the evanescent part of MTF with rapid phase modulation 
is filtered out, and the PSF becomes narrower and gains a regular Gaussian shape. 
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